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Size effects in the nonlinear resistance and flux creep in a virtual 
Berezinskii-Kosterlitz-Thouless state of superconducting fllms. 
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We show that the size effects radically affect the electric field-current {E — I) relation of super- 
conducting films. We calculate E{J) due to thermally-activated hopping of single vortices driven by 
current / across the film in a magnetic field H , taking into account interaction of free vortices with 
their antivortex images and peaks in the Meissner currents at the film edges. Unbinding of virtual 
vortex-antivortex pairs not only mimics the transport uniform BKT behavior, it can dominate the 
observed E(J) and result in the field-dependent ohmic resistance at small I. We show that E{I) 
can be tuned by changing the film geometry and propose experimental tests of this theory. 

PACS numbers: 74.20. De, 74. 20. Hi, 74.60.-w 



The Berezinskii-Kosterletz-Thouless (BKT) transition 
is a 2D universal phase transition due to unbinding of log- 
arithmically interacting topological excitations [l|. The 
concept of the BKT transition first introduced in the con- 
text of vortices in Xy-magnets has been extended to 
other topological excitations like vortex-antivortex pairs 
in superfluid filins, superconducting films, Josephson- 
junction arrays [2i, 3] , dislocation pairs in the theory of 2D 
melting or ultracold atomic gases in optical lattices [J]. 
The superconducting films and Josephson arrays have be- 
come the main experimental testbeds to study the BKT 
transition by dc transport measurements. In this case 
the ohmic electric field-current characteristics R = RI 
above the transition T > Tbkt turns into the power-law 
E oc /!+" at T < Tbkt with a jump to a = 2-5 followed 
by the growth of a as the temperature T decreases [3| . 

While the interaction of dislocations and vortices in 
Xy-magnets of superfluid films is indeed logarithmic, 
the interaction of vortices in superconducting films is 
only logarithmic over distances shorter than the Pearl 
screening length A = A^/d where d is the film thick- 
ness and A is the London penetration depth [5|. The 
size effects can change the BKT transport behavior at 
T < Tbkt since the result, E ex /^+", holds only at 
sufficiently high currents, I > Ii '^ ce/0O: for which 
the critical size of a dissociating vortex-antivortex pair, 
ic = 2cwe/(j)oI, is smaller than the film width, w, where 
a = 2e/T, e — ^q/IGtt^A is the vortex energy scale, (po is 
the fiux quantum, and c is the speed of light. For I < U, 
the E — I characteristic becomes ohmic [2, y, la, [Zl, |8|. 
Yet several crucial features of the electrodynamics of su- 
perconducting films have not been incorporated into the 
BKT theory. First, the sheet current density J{x), which 
drives vortices across the film can be highly nonuniform. 
For a current-carrying thin film strip of width w > A in 
a perpendicular magnetic field H , we have [9[: 




J{x) = [/ + (w - 2x)Hc/A]/T:yJx{w-x), (1) 
where the geometry is shown in Fig. la. This distribu- 



FIG. 1: A thin film in a perpendicular field H . The black dot 
shows a vortex moving across the film, and the empty circle 
shows the antivortex image (a). Geometries for probing the 
resistive state by relaxation measurements: a thin ring in a 
perpendicular field (b), and a thin tube on a cylindrical sub- 
strate in a parallel field (c). The white dot shows the vortex 
driven along the tube by the azimuthal Meissner currents. 



tion of J{x) ensures no spontaneous vortices generated 
by small / and H in the film (the singularities at the 
film edges are cut off at the distances ~ max(d. A)). The 
second feature results from the Bean-Livingston surface 
barrier: a vortex penetrating a film interacts with a ficti- 
tious antivortex image, which provides zero normal cur- 
rents at the edges. Thus, thermally-activated penetration 
of single vortices is governed by the BKT-type unbinding 
of a virtual vortex - antivortex pair [6|. For w < A, the 
interaction energy /7(ri, r2) between two vortices is loga- 
rithmic only for small separation, |ri— r2| < w, otherwise 
^(ri, r2) decays exponentially over the length w/tt along 
the film because of cancelation of the vortex currents by 
an infinite chain of vortex-antivortex images [10| . This 



makes rare therraally-activated hops of vortices across 
the strip uncorrelated at low T and /. 

In this Letter we show that fluxon hopping mediated 
by the unbinding of a vortex from its edge antivortex 
images mimics the uniform BKT resistive state and re- 
suhs in a strongly size-dependent E{I), which can exceed 
E2{I) caused by the uniform pair dissociation [3| both 
for w < K and w > K. This is due to the fact that the 
energy activation barrier for the single vortex penetra- 
tion is roughly half of the barrier required to create a 
vortex-antivortex pair in the film. The account of these 
features is important for the interpretation of deviations 
from the BKT scenario and critical currents observed on 
E - I curves of ultrathin films [O, [HE Q] ■ Since it 
is the thin film strip geometry, which is mostly used in 
dc transport measurements, we also discuss other geome- 
tries in which the genuine BKT pair dissociation could 
be revealed. 

We calculate E{I) due to vortex hopping across a thin 
film described by the Langevin equation rjx -\- U'{x) — C 
where the dot and the prime denote differentiations over 
time t and coordinate x, respectively, rj is the viscous 
drag coefficient, C(i) describes thermal noise and the local 
energy U{x) = Uo—Um comprises the position-dependent 
vortex self-energy Uo{x) and the work of the Meissner 
current, Um = {4'o/c) Jq J{u)du to move the vortex by 
the distance x from the film edge. Here J[x]^ i_s_ described 
by the integral Maxwell-London equation 



x) i 

SI 



J{u)du 



4:TTAd,:J = ~cH 



(2) 



supplemented by the condition I ~ Jn J{x)dx. If w 3> 
A, Eq. ^ yields Eq. ([T|), but for u; <C A, the integral 
term is negligible, and J{x) ~ I /w + cH{w — 2a;)/87rA. 

The self-energy Uo{x) = — /„ F{u)du is the work re- 
quired to create a vortex at the edge where t/o(0) = 
and move it by the distance x. Here F{x) — f{2x) + 
^^j^[/(2K;n -I- 2x) — f{2wn — 2x)] is the force between 
the vortex in the film and an infinite chain of vortex and 
antivortex images outside the film, f{x) — 4>Qjy{x)lc, 
and Jy{x) is the y-component of the sheet current den- 
sity of the Pearl vortex in an infinite film. Using Jy{k) = 
~ic4>akx/2'Kk{l + 2Afc), k'^ = k^ + ky 5] and integrating 
over ky in the Fourier space, we obtain: 



(/>Q Y^ sin {TTnx/w) 
° " ^ ^^ v/(27rnA)2 - w 



: tan 



27rnA 



-,1/2 



27rnA 



(3) 

Here N ^ we *^/27r^ and C ~ 0.577 provide the vortex 
core cutoff. For narrow films w <C 27rA, the summation 
in Eq. ^ reproduces the known result [lO|, ll5( : 



Uo{x) — eln[(z«/7r^) sin(7ra;/u')]. 



(4) 



where e 



(J)q/16tt'^A. Here Uq results from the kinetic 
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FIG. 2: The vortex energy U{x) given by Eqs. ^ and ^ 
for a strip with w <^ A and w = 20^. The dashed line shows 
Uo{x), and the solid lines show U{x) for H(j)ow^ /SivAe — 9 
and different currents, (pol/ce: (1); 1 (2), and 2 (3). 



distance ~ ^ from the edges where the London theory 
breaks down. For wide films w > 27tA, Uo{x) increases 
from zero a.t x = to Ua — eln(A/^) over the length 
a; '--^ A. The magnetic part of the energy barrier Um{x) 
for w 3> A and w ^ A, is given by 



Um = sm ^/x/w + 



Jxiw — x), w > A(5) 
Um = (jiolx/cw + 4>qHx{'w — a;)/87rA, w <^ A(6) 

The behavior of U{x) at different / and H is shown in 
Fig. 2. The transport current tilts U{x), reducing the 
barrier maximum and shifting its position xo{I) toward 
the film edge. The barrier disappears at I ^ Is for which 
xa{Is) ~ £,■ In turn, the magnetic field at / = leaves 
U{x) symmetric, but can produce a minimum in U{x) at 
X — 'w/2. There are 3 characteristic fields: Hi, at which 
the minimum in U{x) appears, the lower critical field 
Hci at which U{w/2) — 0, and i/^, at which the edge 
barrier disappears. These critical currents and fields can 
be calculated from the equation U'{xq) = 0. 

We calculate E(I) for T < Tbkt, H < Hci{I) 
and / < /s(-ff) so that the voltage V results from 
thermally-activated hopping of vortices and antivortices 
over the barrier U^{x) = Uq — U^. Here U^{x) for 
antivortices is given by Eq. ([5]) with H —>■ —H and 
sm~^{x/wy/^ -^ cos~^{x/wy/^ or by Eq. ([6]) with 
H —t —H and x — > w ^ x. The mean drift velocities 
v± of vortices and antivortices follow from the solution 
of the Fokker-Planck equation with a constant probabil- 
ity current [l6[: 



■qv± 



wT\F±{w)-F±{0)]j^ 



"'My) 



dxF±{x) X 

F±{wy 



dy 



F±iy) 



(7) 



energy of unscreened vortex supercurrents cut off at the 



where /3 = e/T, F±{x) = exp[(;7±(a:) - Uo{x))/T], so 
that F+(0) = F^{w) = 1, and F+{w) = F_(0) = 



exp{(f)oI / cT) . The integral over x is cut off on the scales 
of the vortex core size, and the condition T < Tbkt 
implies that /3 > 2. If / ^ J^, where Ig for w <^ K 
is of the order of the depairing current, the x-integral 
is determined by the vicinity of the edges. Indeed, for 
x « 0, the self-energy Uo{x) ~ eln(a:/^) is dominated 
by interaction of the vortex with the nearest image, thus 
F{x) ~ {^/xYF{Q), the first y-integral in the brackets is 
negligible and the lower limit of the second y-integral can 
be set to a; = 0. Doing the same for a; « w, we obtain 
the factor 2F(0)F(w)^/(/3 — 1) after integration over x. 
The velocities v± are proportional to the mean electric 
field E ~ (t>a{v+ — v-)/'wcS,. This follows from the Joule 
power IV = 4>oI{v+ — v-)L/^wc produced by the driving 
force I<J)q/wc to move a vortex across the film and mul- 
tiplied by the number ~ L/^ of statistically-independent 
edge sites available for uncorrelated vortex entries in the 
strip of length L. Using the Bardeen-Stephen expression 
for rj ~ d4)l/2TT(^c^p.n in Eq. (O, we obtain [l7| 



E 



lI2PllIii^}l[l_^-iMcT^[Z^i^ZZ'], (8) 



Z+ = / e^ /^dx. 



(9) 



The behavior of E{I, T, H) described by Eq. ^ is shown 
in Fig. 3: E{I) is ohmic for IcfiQ <^ cT and nonlinear at 
higher /. The ohmic E = R^I at i/ = is quantified by 
the Arrhcnius-type resistance Ry oc (—Ua/T) per unit 
length, for which Eqs. (g]), (g]) and dH) give: 



Ry 



dwT[{l3 - l)/2] 



w < A 



(10) 



where T{x) is the gamma-function. The barrier hight, 
Ua = eln(w/7r^) = Uo{w/2) depends logarithmically on 
w in accordance with Eq. (JH). For /3 ^ 1, Eq. (llOp yields 
Ry ~ {y/2pn / dw){TT P)"^/"^ {-kS, / wY ^ much smaller than the 
normal resistance i?„ = pn/dw. In wide films w S> A, 
the barrier Ua — eln(A/^) becomes independent of w. 

For 2(f)Ql > TTcT, or (p^Hw > 2ttT the change of the 
barrier shape U{x) shown in Fig. 2 results in a strongly 
nonlinear and field dependent E{I, H), which can be cal- 
culated numerically from Eqs. ([2]) and ([8]) for any ratio 
w/A, and analytically for both limits w ^ A and w ^ A. 
For instance, in wide films at Is£,/\/dw <C / <C Is, the 
fiuxon hopping is limited by the small barriers near the 
edges: U+{x) ~ eln{x /£_)- (f>o{2I / nc + IIw/2n)y/x/w a.t 
X <^ w. For H 3> 2'kT/w(I)o, the antivortex channel is 
suppressed, Z+ <C ^-, so Eqs. ^ and Q yield: 



E 



7:pncT{p~l){^/wy 
2dw(f)oTi2P + 2) 



^/2/ 
ttT \ c 



Hw 



213+2 



(11) 

In the limit /0o < cT, but II4>qw > 2tiT, the ohmic 



resistance Ry strongly depends on H: 



Ry 



TTPnjP - 1) 

2dwr{2/3 + 2) 



hHw\ 
2ttT ) 



2,9+2 



(12) 



For H ^ 4:1 /c, but 2I(f>Q 3> ncT, the vortex and antivor- 
tex channels yield the power-law E{I): 



E ^ 



np^cT{(i-l) (^y (2<j>oI 



dw4)QT{2(3 + 2) \wj \ncT 



2/3+2 



(13) 



For narrow film w ^ A at H — 0, the integral in Eq. Q 
can be evaluated analytically for all I < Ig-. 



E = 



AirpncTiP - 1) 
d(l)owT{l3 + 1) 



27rf 



2 

r(l-f-^-fi7)| sinh7r7 (14) 



where 7 = 4iqI/2ttcT. In the limit 7 <C 1, Eq. ^ 
reproduces Eq. PU)) . but for 7 » /3/2, that is, JqS./w «; 
J < Jf) where Jq = ccfio/Sn'^eA^^ is of the order of the 
sheet depairing current density, Eq. (|14p gives 



E = 



27rp„if3-l) f MJ 
dT{(3 + l) \ cT 



J. 



(15) 



This power-law E{J) can also be obtained in the same 
way as Eq. (|13p by expanding U{x) near the film edges. 
Notice that E{J) given by Eq. (fT5|) is independent of w 
because, once the vortex overcomes a narrow (<C w) edge 
barrier shown in Fig. 2, its subsequent viscous motion 
across the film is no longer thermally-activated. 

It is instructive to compare Eqs. p3|) and p5|) with the 
electric field £'2 ^ (p„ J/d)(J/Jo)^'^ produced by the uni- 
form BKT dissociation of vortex-antivortex pairs above 
the critical size £c — 2ec/(j)oJ |2]. For narrow films 
at low temperatures, (w ^ A, /3 » 1), we can use 
r(z) ~ (27r/z)i/2e-^z^ in Eq. ^ and obtain 



E2/E ^ {J/2Jof/ey 



(16) 



Hence, for /? 3> 1, the virtual vortex-image unbinding 
dominates over the uniform pair dissociation except in 
the region T w Tbtk of the genuine BKT behavior. 
In wide films, the single- vortex contribution E/E2 ~ 
i^/O^ ^ 1 is further enhanced by the singularities of 
the Meissner current at the edges. As an illustration 
Fig. 3 shows E{I) calculated from Eq. ([8]), which gives 
E > (10^ — 10'^)i?2 in the region where ic < w. Moreover, 
E{I) due to the edge vortex-image unbinding exhibits all 
characteristic features of the BKT nonlinear transport in 
a finite size film: the ohmic E{J) below the critical cur- 
rent Ic followed by the power-law E = RJ{J/Jo)°'^ for 
I > Ic- Here the exponent ai varies from ai =2(3 + 1 
for wide films to ai = /3 for narrow films, while the uni- 
form pair unbinding gives 02 = 2/3 [2]- The similarity of 
ai and a2 in wide films results from the edge Meissner 
singularity of J{x), which increases ai as compared to 
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FIG. ?,: E - 3 curves calculated from Eq. © for /3 = 4, 
It) = 20^ < A and different fields h = J/^o^VSttAT. 
Here j = /(jio/cT < js = 2f3w/e£_, and Eq = TvcpnT{l3 - 
1){tv^/w)^ /dw(f>o. The line labeled by HN shows the Halperin- 



Nelson result, E2{I) Q]- The critical pair length . 
in the region j < 2/3 left of the dashed line. 



exceeds w 



ai = (3 for a uniform J . The critical current Ic is esti- 
mated from the condition that the maximum of U{x) at 
X = xo(/) shifts from the fihn center at I <^ Ic to the 
edge at xo ^ w for I ^ Ic- For a narrow film, Ic defined 
by xq{Ic) = w/4 in Eqs. ^ and ^ is: 



Ic{H) 



IGttA 



1- 



i/ 
i^ 



Ho = 



(17) 



so that /c(0) is independent of w, but both Jc(0) = 
/c/w ~ Jq^/ui and -ffg increase as w decreases. The same 
/c(0) = TTce/fjJo is obtained, defining the nonhnearity on- 
set from the condition Ic = 2w/n equivalent to 7 = /3/2 
in the argument of the gamma- function in Eq. (|14p . 

The results presented above indicate that E{J) can 
be tuned by changing the film geometry. For instance, 
if a uniform J{x) is produced in a wide film, the expo- 
nent ai =2(3+1 would decrease to ai = /3. This could 
be implemented by using ferromagnetic/superconducting 
structures 15|, [l8[, in which a thin film strip is placed 
perpendicular to ferromagnetic screens to eliminate the 
singularity in J{x) 18]. Another possibility is to use a 
thin film tube in a parallel field, which produces uni- 
form azimuthal screening currents J = cHd/AirX driving 
vortices along the tube. Because of the negligible demag- 
netization factor of this geometry, J{x) for large tubes of 
length L ^ A and diameter Z? 3> A does not contain the 
Meissner edge singularities characteristic of wide films in 
a perpendicular field. Such a tube would have a mixed 
E{J) controlled by Uo{x) of a wide strip, but a uniform 
current drive like in a narrow film. 

Film and ring structures make it possible to probe 
E(J) by magnetic relaxation measurements well below 



the nV voltage sensitivity ll| of transport measure- 
ments. In this case i?(i) is ramped up and then stopped, 
after which the magnetic moment M{t) = I{t)D/2c 
is measured. For £1 3> 0oc, relaxation of I{t) in 
a ring or a tube is described by the circuit equation 
£1 = —TTc^DRI{I/Io)", where £ is the self-inductance. 
The solution of this equation, I{t) = {T/t)°'Io with 
r = £/Trc^DRa, enables extracting a{T) from fiux creep 
measurements after some initial transient time [19|. 

In conclusion, thermally-activated fiuxon hopping me- 
diated by unbinding of single vortices from their edge 
antivortex images can mimic the nonlinear resistive be- 
havior of a uniform BKT state. Our results predict a 
strong dependence of E{J, H, T) on temperature, mag- 
netic field and the sample size. This offers a possibility 
of tuning the behavior of E{J) by changing the film ge- 
ometry or by incorporating magnetic structures. 
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